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Abstract 

We examine the possibility of determining the neutrino mass hierarchy in the limit 
013 = using atmospheric neutrinos as the source. In this limit, in which ^13 driven 
matter effects are absent, independent measurements of A31 and A32 can, in principle, 
lead to hierarchy determination. Since the difference between these two is A21, one 
needs an experimental arrangement where A21L/E > 1 can be achieved. This condition 
can be satisfied by atmospheric neutrinos since they have a large range of energies and 
baselines. In spite of this, we find that hierarchy determination in the ^13 = limit 
with atmospheric neutrinos is not a realistic possibility, even in conjunction with an 
apparently synergistic beam experiment like T2K or NOi^A. We discuss the reasons for 
this, and also in the process clarify the conditions that must be satisfied in general for 
hierarchy determination if ^13 = 0. 
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1 Introduction 



The neutrino mass hierarchy is an important discriminator between large classes of models of 
unified theories [1]. Its determination remains one of the outstanding problems facing neutrino 
physics. The term hierarchy refers to the position of the third mass eigenstate with respect 
to the other two (relatively) degenerate states. When the hierarchy is normal (NH), the third 
mass eigenvalue {m^) is the largest and hence A31 = mg — ml and A32 = mg — m"^ are both 
positive. On the other hand, when it is inverted (IH), m-^ is the smallest, leading to A31 and 
A 32 being negative. 

The most widely discussed method for hierarchy determination relies on sizable matter 
effects at long baselines. The baselines of the foreseeable super-beam experiments are 1000 Km 
or less. Therefore in these experiments appreciable matter effects do not develop, although it 
may be possible to have some hierarchy sensitivity in the z/^ — > i/g channel [2,3] provided ^13 
is large enough (z.e.sin^ 2^13 > 0.04). However, this channel is seriously compromised by the 
((5cp, hierarchy) degeneracy [4,5]. Thus hierarchy determination via this channel requires data 
from a number of complementary experiments. The degeneracy problem is reduced in the 
case of a detector located at the magic baseline of ~ 7000 km [5-7]. However, this requires 
high luminosity sources such as neutrino factories or beta-beams [8], likely to be available 
only in the far future. Recently synergy in T-conjugated channels were explored to determine 
hierarchy for shorter baselines with small matter effects for sin^ 26'i3 ^ 0.03 [9]. 

Independent of the value of 6*13, A31 is the largest mass squared difference for NH, whereas 
for IH, A32 = A31 — A21 has the largest magnitude. In principle, one can exploit this difference 
to determine the mass hierarchy if one can measure | A31 1 and | A32 1 individually and precisely 
[10]. Detection of the interference effect in vacuum neutrino oscillations [11-13] and the 
measurement of the phase of monoenergetic [14] have been proposed to determine the 
hierarchy. 

With atmospheric neutrinos as the source, wide ranges in energy (E) and baseline (L) 
become available and it is possible to observe large resonant matter effects [15-20]. For values 
of 6*13 much below the CHOOZ upper bound [21], the distance a neutrino needs to travel for 
these effects to be observable, becomes larger than the diameter of the earth [16]. There- 
fore, hierarchy determination with atmospheric neutrinos using matter effects also requires a 
moderately large value of 6113 [18,19] (sin^ 26'i3 ^ 0.05). 

There is a large class of models which predict 613 to be exactly zero due to some symmetry 
principle [22] 0. Recent global analysis of neutrino data indicates a small preference for a non- 
zero 6*13, sin^ 26*13 — 0.05 [30,31]. However, at present the statistical significance of this result 
is only at the level of ~ la. In addition, the latest results from MINOS show an excess of 
electron events [32]. This may also be an indication for a non-zero ^13 [33]. 

In the present work we consider the problem of determination of the sign of atmospheric 

^It is of course likely that such symmetres are broken at low energies and ^13 would acquire a small non- 
zero value [23-27]. Hierarchy determination using matter effects driven by very small values of ^13 entails 
measuring the Vf, probability at the magic baseline using a beta-beam for sin^ 2(?i3 > 10^^ [28], or with 

a neutrino factory beam for sin^ 2^13 > 10""^ [29]. 
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neutrino mass difference in the limit ^13 — > 0. This is a challenging problem requiring high 
precision measurements. The general requirements for this were discussed in [34]. It was 
argued that one would need measurements of P^^ and P^p at different L/E with atleast 
one L/E satisfying A21L/-E' > 1. In this conenction the utility of a broadband accelerator 
beam or atmospheric neutrinos capable of providing different values oi L/E within the same 
experimental set up was pointed out in [34]. Reference [35] made a detailed implementation 
of the ideas discussed in [34] in the context of neutrino factories and superbeams using two 
narrow band beams at different L or one broad band beam. However, no detailed quantitative 
analysis of hierarchy determination for ^13 = using atmospheric neutrinos had been done 
in the literature so far. In this paper, we undertake this task and perform the analysis 
using atmospheric neutrinos in conjunction with long baseline experiments. Our study shows 
that despite the fact that atmospheric neutrinos provide a large range of L and E, practically 
speaking, this is an impossible prospect. We discuss in detail the conceptual and experimental 
issues that lead to this conclusion. 



2 P^^ in the limit ^13 = 

For ^13 = 0, the probabilities involving the electron flavour depend only on A21 = — m\ 
because the mass eigenstate decouples from z/g. Hence probabilities involving Ug. are not 
sensitive to the sign of A31 and to determine the hierarchy, other oscillation and survival modes 
must be explored. In the case of atmospheric neutrinos, the z/^ — >■ v^^ survival probability 
provides the most effective channel, especially from the detection viewpoint. 

As explained above, for different hierarchies the relative magnitudes of A31 and A32 are 
different. The difference in magnitudes of A31 and A32 for NH and IH is proportional to A21. 
To exploit the corresponding difference in P^^, one needs values of L/E such that A21///-E' > 1 
or L/E > 10^. Such values of L/E can be obtained for E <1 GeV and L ~ 10^ Km. A large 
flux of atmospheric neutrinos with these baselines and energies is readily available. 

For such long distances, the neutrinos propagating through earth's matter experience a 
potential V = \plGFne [Gp is the Fermi Constant and is the ambient electron number 
density). This gives rise to an effective mass squared matter term, A = 2EV which can be 
expressed as 

A = 2 X 7.6 X 10"^ p(g/cm3) E{GeV) [eV^] (1) 

where p is the average density of matter along the neutrino path and Yg = 0.5 is the electron 
fraction per nucleon. 

Including matter effects, the expression for P^^ for ^13 = is 

P- = l-4cos^^23sin2r2Cos2r2sin^(Ari^) 

-4 cos^ ^23 sin^ ^23 sin2 9^2 sin2 A™ 

-4 cos^ 023 sin^ 023 cos^ sin^ (^^2^) (2) 
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where, the squared mass-differences in matter are 
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v/( A21 cos 2^12 - Ay + ( A21 sin 2^12) 

sin2er2 = sin 2^12^. 

'^21 

From the above equations we see that the condition A21 cos 26*12 
energy related to the solar mass-squared difference A21, 

^i?/^ A21 COS 26*12 



(3) 

(4) 
(5) 

A defines a resonance 



(6) 



0.76 X 10-^ p{g/cm^)' 

The present best fit value of A21 ~ 8 x IQ-^ eV^ [36-38] gives ^ 0.06 - 0.2 GeV. 

In the third column of Table 1 we list the values of for the various baselines. Note 
that for longer baselines, the average density is larger because the neutrinos pass through the 
inner mantle (5000 < L< 10000 Km) and the core (L > 10000 Km). 



L {km) 


Pavg {gm/cm^) 


E^ (GeV) 


295, 732 


2.3 


0.20 


2900 


3.3 


0.14 


7330 


4.2 


0.11 


12000 


7.6 


0.06 



Table 1: Values of arc listed as a function of baseline or the density p. See text for details. The value of 
sin^ 2^12 = 0.8 or cos 2^12 = 0.43 is used for evaluating E^. 



For energies E « E^, the matter term is negligible and we obtain the vacuum limit of 
the survival probability. Its form is similar to eq. ([2]) with the vacuum angles and mass-square 
differences taking the place of the corresponding matter dependent quantities. For energies 
E » or A » A21, we get 



A^} ^ A - A21 cos 2^12 and ~ 7r/2. 
Substituting these in eq. ([2]), we get 

L 



(7) 



1 — sin 2^23 sin 



(A31 - A2l4) 



AE 



Here we define c,- 



cos 6ij and Sij 



sin 6ij. We note that the survival probability in matter is 



independent of the matter term because of cancellations [34] . But it must be stressed that this 
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probability is in general not the same as the vacuum survival probability. It coincides with 
vacuum P^^ only in the limit A2i-L/4£' <^ 1 and if only terms linear in the small parameter 
A21/A31 are retained in both expressions. In this limit both P^^ in vacuum and P™^ from 
eq. (IH]) are given by 



J- nn — 



1 - sin^ 2^23 



sm 



A31L 
AE 



AE 



sm 



A31L 
2E 



(9) 



It is to be noted that the lengths involved in the long baseline experiments T2K [39] (L = 295 
Km) and MINOS [40,41] and NOi^A [42,43] (L = 732/810 Km) are short enough such that 
the approximation A21L/4E <^ 1 is valid. The energy of T2K {E > 0.5 GeV) is a few times 
E^ whereas the energies in MINOS and NOz/A [E > 2 GeV) are much larger than E^ as 
can be seen from the first row of Table 1. Therefore for MINOS and NOi^A eq. ([H]) is valid 
to a very good approximation. It is interesting to consider the precise difference between 
measured by T2K [39] and the corresponding vacuum probability in eq. ([9]) . Expanding terms 
in eq. ([2]) to first order in the small parameter A21L/4E, we get 



P. 



1 - sin^ 29. 



23 



sm 



Am T 
AE 



{c^2?A-[L ^.^ A-L 



AE 



2E 



(10) 



In this expression, (c^)^ is very small (~ 0.09) because 6^ is just a little lower than 7r/2. This 
multiplies the small term A^^L/E (which is of the same order as A21L/E). Hence the second 
term in the square bracket can be neglected and we have 



pm 



1 



A"^L 

sin^ 2^23 sin^ 



1 - sin^ 26'23 



sm 



AE 
A31L 
AE 



O.55A21L . A31L 
sm 



AE 



2E 



Comparing this expression with eq. above, we find that the only difference is the replace- 
ment of the factor c^2 — 0-69 by 0.55 in the coefficient of the second term. But since this 
term is suppressed by the small parameter A21L/AE = 0.04, compared to which the contri- 
bution of the first term (sin^ A3iL/4£') is about 0.9 for T2K, the difference in the coefficient 
causes a negligible change in the value of P^^. The difference in magnitude between P^^ from 
eq. and eq. (ITTI) above is less than 0.5%. This makes eq. (ITT]) equivalent (to a very good 
approximation) to eq. (E]) in this limit [34]. 



3 Hierarchy Sensitivity 

It is clear from eq. ([H]) that if the magnitude of (A31 — A2icf2) is different for NH and IH, P^^ 
will also be different for the two hierarchies. Before computing the difference in P^^, we need 
to figure out the difference in the magnitude of the quantity (A31 — A2icf2) with a change 
in the hierarchy. We investigate this change in the light of various assumptions made in the 
literature. 
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• First we consider A31 (Jif) = —A^i^NH). This assumption amounts to identifying A31 
with Am^^^ determined by Super-K [44] and is widely used in the hterature for situations 
where the difference between IA31I and IA32I hes well below experimental errors, and is thus 
immaterial. It leads to 

p;;(NH)-p-(iH) = 

sin2 2^23 [sin^ [(A31 + A2iC?2)L/4E] 

- sin2[(A3i-A2i4)L/4E]]. (12) 

Clearly, this assumption is untenable for the situation under consideration here, where 
the difference between these two mass differences is crucial. 

• The other assumption often used is A3i(/if ) = -A32 (iVif) = - A3i(iViJ) + A21. This is 
equivalent to the statement that the largest mass squared difference has the same magnitude 
for both NH and IH. It leads to 

(A3i(/if) - A2i4) = {-Au{NH) + A2i(l - 4)) 

= -{Asi{NH)-A2isl,). (13) 

Substituting this in eq. ([H]) gives 

p;rM(NH) - p;;(iH) = 

8111^2023 [sin^ [(A31 - A2i4)L/4E] 

- sin2[(A3i-A2i4)V4^]]- (14) 

This amounts to making an ad hoc assumption which currently is unsupported by exper- 
imental evidence. 

• If we do not make any assumptions but instead consider the question: what information 
do the experiments give regarding the magnitudes of A31 and A32 for NH and IH ? Ongoing 
experiments such as MINOS and future experiments such as T2K and NOi^A measure the 
muon neutrino survival probability. In the two flavour limit, 613 = and A21 = 0, this is 
given by 

P,V' = l-sin^2^23sin^(^). (15) 

Analyzing the data of these experiments, in this limit, gives us the magnitude of A31. This 
is the reason behind the widely used assumption A31 (/_?/) = —A^ii^NH), quoted at the 
beginning of this section. 

The survival probability for T2K and MINOS/NOi^A is given by eq. (IH]) to a very good 
approximation as shown in the previous section. The three flavour effects appear in the form 
of non-zero A21. We note that eq. ([8]) can be obtained from eq. (fT5|) by the replacement 
A31 — i> (A31 — c^2^2i)- Therefore the "atmospheric" mass-squared difference measured by 
MINOS or T2K is not A31 (or A32) but the combination [13,34]. 

AmL, = A31 - 4A21 = sl^Asi + cl,A32. (16) 
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Thus, in the hmit = 0, MINOS or T2K measure the magnitude of the mass-squared 
difference Am^j„^ defined in eq. fll6p . This quantity is positive for NH and negative for IH. 
Then from eq. (fT6l) one can easily derive A^ilIH) = —A^i{NH) +2c^2^2i- With this relation 
between A^ilNH) and A^i{IH), we see that P™^ in eq. (IHl) has no hierarchy sensitivity [34]. 
Atmospheric neutrino data with energies E » Ef measure the above P™^, and hence contain 
no hierarchy sensitivity. This statement is true independent of the precision with which Am\^^ 
can be measured. Any deviation from this prediction is likely to come from the regions where 
the approximations made in obtaining eq. (E]) are not exact [34,35]. This is equivalent to saying 
that a combination of two experiments can give hierarchy sensitivity only if the mass-squared 
differences (i.e. the frequencies) measured by the two experiments are different. 

In Figure [H the hierarchy difference for P™^ with 6*13 = is plotted as a function of the 
neutrino energy E for a baseline of 12000 Km using the three definitions of the hierarchy 
discussed in this section, i.e. (a) A31 (JiJ) = -A-ii{NH), (b) A31 (JiJ) = -A-^^NH) + 
A21 and (c) A31 (Ji7) = —A3i{NH) + 2c^2^2i- Note that the assumptions a) and b) lead 
to a large but spurious sensitivity to the hierarchy, as reflected in the large values of the 
hierarchy difference for P™^. This difference is greater for ) than case (b). Case (c) 

gives the hierarchy difference for P™^ for the third hierarchy definition, which is based on 
the experimental measurement of Am^^^ by accelerator experiments. This implies that we 
need this measurement from experiments such as T2K or NOi^A to explore the hierarchy 
sensitivity in atmospheric neutrino data. In generating the plots, P/^/i is obtained by the 
numerical integration of the neutrino evolution equations. 

We see that the hierarchy sensitivity is quite small for case (c). A difference ~ 0.2 can be 
observed in the numerical plots at lower energies (i? ~ 1 GeV) for this very long baseline of 
12000 Km. This is because the expressions given in eq. and eq. ([H]), derived in the constant 
density approximation, do not hold good at such values of E and L. But this difference is very 
small and would be washed out even with optimistic values of energy and angular resolution 
of the detector. 

From Table 1, we see that for multi-GeV atmospheric neutrinos the energy E will always 
be larger than E^ and so it will be difficult to find a regime where the probability will be 
different than that given by eq. Table 1 also shows that E^ is larger for shorter baselines. 
Thus, we expect that a departure from the form given in eq. (IHl) is most likely for the baselines 
of 3000 — 6000 Km and neutrino beam energy of 0.2 — 0.7 GeV. These scenarios were analysed 
in [35] for superbeam and neutrino factory experiments. 

If we look at the region in atmospheric neutrino experiments where the neutrino energy is 
moderately larger than E^ (in the range 0.2 - 0.6 GeV), the matter modified mixing angle 6^ 
will be less than 7i/2. Hence, we can't use P^^ in eq. but need to go to the full expression 
given in eq. ([2]). This equation is in principle sensitive to hierarchy due to the matter term A 
and the difference between A31 and A32. However, this difference oscillates rapidly with energy 
and with baseline. To measure it, one would require extraordinary neutrino energy and angular 
resolution and very large statistics. Such angular resolution is impossible to achieve in the 
relevant sub-GeV energy range, because quasi-elastic neutrino scattering, which has a broad 
distribution in the direction of the final state lepton, dominates in this range. Additionally, 
most atmospheric detectors with target nuclei have an inherent limit on the energy resolution 
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possible, set by Fermi motion of the bound nucleons. Such energies are typically ~ 100 MeV, 
and thus are a significant fraction of the energy of the final state lepton in the range of interest 
here. 

In Figure 2, we plot the difference in P™^ for NH and IH as a function of baseline for the 
two energies, 0.5 GeV and 1.0 GeV. These figures show a moderate to substantial difference 
for the case of no smearing. However, this difference becomes less than 0.02 when we include 
an angular smearing with = 5°. This value of ag is much smaller than the width of the 
angular distribution for quasi-elastic events. Note that no energy smearing was included in 
generating this figure. Thus we see that even with an ideal energy resolution, the difference 
in P™^ for the two hierarchies becomes tiny. If a moderate energy smearing is included the 
difference is likely to vanish. Hence it is not possible to determine the neutrino mass hierarchy 
using atmospheric neutrinos as a source, no matter how good the detector is. 

To summarize, the general requirements for hierarchy determination for ^13 = are 

1. |Am^j^| should be measured to a precision better than 2%. 

2. The neutrino energy E should be large enough to produce a muon but should not be 
too large compared to E^. 

3. The baseline L should be such that A21L/E > 1. 

4. Excellent energy and angular resolution are necessary, so that rapid oscillations in P^^ 
can be resolved. 

Atmospheric neutrinos clearly do not satisfy the last criterion. In a beam experiment, 
however, ae is zero, so the problem of angular resolution is automatically resolved. However, 
it must still meet the requirement of superior energy resolution, good enough to resolve the 
closely spaced peaks shown in case (c) of Figure 1. It is apparent from the size of the peaks 
that the experiment needs to have the capability to gather high statistics [35]. 

4 Conclusions 

In conclusion, we have studied the feasibility of determining the sign of A31 for 613 = using 
atmospheric neutrinos as the source. When Oi^ = 0, matter effects induced by it are absent. 
Moreover, z/g is decoupled from u^, leading to no hierarchy dependence in Pee or P^e- This 
makes P^^ the most suitable remaining channel for hierarchy determination. In this limit, 
depends on A21, A31 and A32. For NH (IH), A31 (A32) is the highest frequency. If these 
two close frequencies can be resolved along with their magnitudes, then the hierarchy can be 
determined. This requirement is, in practice, difficult to achieve by measurement at a single 
L/E. Hence it translates to having two experiments, one of which measures a single frequency 
(I Am^^^l), which is a known combination of the two independent frequencies (A21 and A31), to 
high precision. The other experiment then needs to satisfy the requirement A21L/E > 1 and 
measure a different combination of these frequencies, also with high precision. The precision 
requirements imply that the energy and angular resolutions have to be exceptionally good. 
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The A21L/ E > 1 requirement implies long baselines. The condition for measuring a different 
frequency enforces tapping into A21 driven matter effects which in turn needs energies which 
are not too far removed from the solar resonance energy Ef. 

Atmospheric neutrinos offer a broad range oi L/E values including ^2iL/E > 1. But 
the precision requirements preclude any chance of hierarchy determination by a single at- 
mospheric experiment on its own. We have explored if hierarchy determination is possi- 
ble using atmospheric neutrinos in conjunction with upcoming accelerator experiments such 
as T2K and NOi/A. For these experiments, L/E values are such that A^iL/E ~ 1 and 
hence A21L/E << 1. In this approximation, P™^ depends on the single effective frequency 
A31 — c^2^2i- For atmospheric neutrinos with energies above 1 GeV and satisfying the con- 
dition A21L/E > 1, one will measure this same frequency. Thus, even if one had a high 
precision atmospheric neutrino experiment detecting multi-GeV neutrinos, this combination 
would not satisfy the requirement of different frequency measurement discussed above. 

If, on the other hand, one had an atmospheric experiment capable of measuring lower 
energies (which would have the consequence of widening the range of L required to satisfy 
A21L/E > 1 and allowing the inclusion of shorter baselines), it is likely that the different 
frequency condition could be met. However, the precision requirements on energy and angular 
resolution become very demanding in this case. In particular, the energy resolution has to 
be better than the typical Fermi energy of a bound nucleon. Also, given that quasi-elastic 
scattering provides the dominant contribution to the event rates for sub-GeV neutrinos, the 
precision requirement on angular resolution is impossible to realize. Hence neutrino hierarchy 
determination using atmospheric neutrinos, even in conjunction with another high precision 
beam experiment, is not possible. 
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